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O . Analysis of the isoscalar S- and D-waves of processes tttt tttt, KK, rjrj 

^ ■ is carried out aimed at studying the status and QCD nature of scalar and 

C^i tensor mesons below 2 GeV and 2.3 GeV, respectively. Assignment of these 

O ! mesons to lower scalar and tensor nonets is proposed. 
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1 Motivation 



The study of spectrum of low-lying hadrons and their properties is very 
important for investigation of the confinement problem and elaboration of 
nonperturbative methods of QCD. Especially it concerns scalar mesons. 
E.g., two different nonperturbative methods of QCD (QCD sum rules and 
unquenched lattice calculations) give essentially distinct results for the 
lightest scalar glueball. 

QCD sum rules: g a lightest scalar meson with a mass below 900 MeV is 
rather narrow and non-qq state (glueball) - see, e.g., [1, 2, 3]. 
Lattice simulations: The lowest mass state of a pure glue should be the 
0++ with a mass of 1670±20 MeV [4]. 

An assignment of the discovered scalar mesonic states to quark-model 
configurations is problematic up to now. 

It is very important to have model-independent information on investi- 
gated states and on their QCD nature. It can be obtained only on the basis 
of the first principles (analyticity and unitarity) directly applied to analysis 
of experimental data. That approach permits us to introduce no theoret- 
ical prejudice into extracted parameters of resonances. We have already 
proposed such method [5] . Here we have applied it to combined analysis of 
experimental data on the processes tttt — ^ tttt, KK, rjrj in the channel with 
the quantum numbers I^J^^ = 0''"0'^"'". In the considered 3-channel case, 
it is turned out to be possible to use a method of the uniformizing variable 
which takes into account the Riemann surface structure. Considering the 
obtained disposition of resonance poles on the Riemann surface, obtained 
coupling constants with channels and resonance masses, we draw definite 
conclusions about nature of the investigated states. 

Further we analyze the same processes in the channel with the quantum 
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numbers I'^J^'^' = 0+2++ for the study of the /2 mesons below 2.3 GeV. 
In this sector, from thirteen discussed resonances, the nine ones (/2(1430), 
/2(1565), /2(1640), /2(1810), /2(1910), /2(2000), /2(2020), /2(2150), /2(2220)) 
must be confirmed in various experiments and analyses [6]. Recently in 
the combined analysis of pp tttt, 7777, 7777', five resonances - /2(1920), 
/2(2000), /2(2020), /2(2240) and /2(2300) - have been obtained, one of 
which (/2(2000)) is a candidate for the glueball [7]. 

In the tensor sector, in addition to the indicated three channels, we 
consider explicitly also the channel (27r)(27r). In the 4-channel case it is 
impossible to apply the uniformizing-variable method with using a simple 
variable. Therefore, the resonance poles are generated by the some 4- 
channel Breit-Wigner forms with taking into account a Blatt-Weisskopf 
barrier factor [8] conditioned by the resonance spins. 

2 Three-coupled-channel formalism 

The S'-matrix for the 3-channel case is determined on the 8-sheeted Rie- 
mann surface. The elements Saf3, where a,/3 = l(7r7r), 2(A'A'), 3(7777), have 
the right-hand cuts along the real axis of the s complex plane, starting 
with 4777,^, 477T,|^, and 47n^, and the left-hand cuts. The sheets of surface 
are numbered according to the signs of analytic continuations of the chan- 
nel momenta 

h = (s/4 - miy/^ k2 = is/4 - 7n|) V2, k^ = (s/4 - Tn^) 
as follows: 
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The resonance representations on the Riemann surface are obtained with 
the help of formulas from Ref. [5], expressing analytic continuations of 
the matrix elements to unphysical sheets in terms of those on sheet I that 
have only zeros (beyond the real axis) corresponding to resonances, at least, 
around the physical region. Then, starting from resonance zeros on sheet I, 
we obtain 7 types of resonances corresponding to 7 possible situations when 
there are resonance zeros on sheet I only in (a) ^n; (b) S22] (c) 533; (d) 
and 5*22; (e) S22 and 5*33; (f) 5ii and 533; and (g) .Sii, S22; and 5*33. E.g., the 
arrangement of poles corresponding to a resonance of type (g) is: each sheet 
II, IV, and VIII contains a pair of conjugate poles at the points that are 
zeros on sheet I; each sheet III, V, and VII contains two pairs of conjugate 
poles; and sheet VI contains three pairs of poles. A resonance of every 
type is represented by a pair of complex-conjugate clusters (of poles and 
zeros on the Riemann surface) . Representation of multichannel resonances 
by the pole clusters gives a main model-independent effect of resonances. 
The cluster kind is related to nature of the state. The resonance coupled 
relatively more strongly to the tttt channel than to the KK and rjr] ones 
is described by the cluster of type (a); the resonance with dominant ss 
component, by the cluster of type (e); the glueball, by the (g) cluster. 
Note that at usual representation of multichannel resonances by the simple 
Breit-Wigner forms, cases (d),(e), (f) and (g) practically are lost. 

We can distinguish, in a model-independent way, a bound state of 
colourless particles {e.g., KK molecule) and a qq bound state [5, 9]. 
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For the combined analysis of data we use the Le Couteur-Newton rela- 
tions [10]: 

^ _ d{-ki,k2,h) ^ _ d{ki, -^2,^3) ^ _ d{ki,k2,-ks) 
d{ki,k2,h) ' d{ki,k2,k3) ' d{ki,k2,k3) 

e c c2 d{-ki,-k2,k3) 2 d{-ki,k2,-k3) 

d{ki,k2,ks) d{ki,k2,ks) 

The Jost matrix determinant d{ki,k2,k3) is the real analytic function 
with the only square-root branch-points at /Cj = 0. 

In the model-independent approach that is based only on the first prin- 
ciples (analyticity-microcausality and unitarity) and is free from dynamical 
assumptions, we use the mathematical fact that a local behaviour of an- 
alytic functions determined on the Riemann surface is governed by the 
nearest singularities on all corresponding sheets. To take into account the 
branch points, we must find proper uniformizing variable. However, it is 
impossible to map the 8-sheeted Riemann surface onto a plane with the 
help of a simple function. Therefore, we neglect the influence of the tttt 
threshold (however, unitarity on the tttt cut is taken into account). This 
approximation means the consideration of the semi-sheets of the Riemann 
surface nearest to the physical region. The uniformizing variable can be 
chosen as [5] 

^2 + k3 

w = I . (2) 

It maps our model of the 8-sheeted Riemann surface onto the w;-plane. 

On the w;-plane, the Le Couteur-Newton relations are somewhat modi- 
fied: 

d*{-w*) d{-w^^) d{w^^' 



Sn — — J]—: — , 5'22 — — — , 5'33 — ,/ X , 
a[W) a[w) a[w) 

c2 _ d*{w*-') 2 _ d*{-w*-') 
d{w) d{w) 



The (i-function is taken as 

d = dsdres 

where the resonance part is 

dresiw) = n + ^r) (4) 

r=l 

with M the number of resonance zeros, ds, describing the background, is 

3 1^ 

ds = exp[-z — («n + iPn)], (5) 

n=l l^n 

where 

Ctn = dnl + ana- —0{s - S^) + ttny- —6{s - Sy), (6) 

S(j Sy 

s — s s — s 

f^n = Kl + Ka -Ois - Sa) + bnv -0(3 - Sy). (7) 

Sfj Sy 

The second terms in a^, and (3n take into account effectively possible chan- 
nels below roughly 1400 MeV (mainly crcr-channel) ; the third terms, many 
opening channels in the range of about 1.5 GeV {r]r]' , pp,ujuj)] Sy is their 
combined threshold. Moreover, the tttt background is taken to be elastic 
up to the KK threshold. 



3 Model-independent analysis of isoscalar-scalar sec- 
tor 

We analyzed the data on three processes tttt tttt, KK, r]r] in the channel 
with the vacuum quantum numbers. As the data, we use the results of 
phase analyses which are given for phase shifts of the amplitudes 6ab and 
for moduli of the S'-matrix elements rjab = \Sab\ (a,b=l-7r7r, 2-KK, S-rjrj): 

Saa = r]aae^''''% = (8) 
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If below the 7777-threshold there is the 2-channel unitarity, then the relations 

iln = il22j ^2 = (1 - ^12 = ^11 + ^22 (9) 

are fulfilled in this energy region. 

For the tttt scattering, the data from the threshold to 1.89 GeV are taken 
from Ref. [11]; below 1 GeV, from many works [12]. 

For TTTT — ^ KK, practically all the accessible data are used [13]. For the 
process tttt — ^ r]!], here we exploited data for the quantity \Sis\'^ from the 
threshold to 1.72 GeV [14]. 

We considered the case with all five resonances discussed below 1.9 GeV. 
From a variety of variants of the resonance representations by possible 
pole-clusters, the analysis selects the following one to be most relevant - 
when the /o(600) is described by the cluster of type (a); /o(1370), type 
(c); /o(1500), type (g); /o(1710), type (b); the /o(980) is represented only 
by the pole on sheet II and shifted pole on sheet III. Description of the 
resonances of types (a), (b) and (c) can be related to the Breit-Wigner 
forms. To reduce the number of adjusted parameters, we make it here, 
except for the /o(980)). 

We obtain a satisfactory description: for the TTTr-scattering from about 
0.4 GeV to 1.89 GeV (xVndf = 202.111/(165-34) ^ 1.54); for the process 
TTTT KK, from the threshold to about 1.6 GeV (xV^df = 161.912/(120- 
33) « 1.86); for the \Si^\'^ data of the reaction tttt 7777, from the threshold 
to 1.72 GeV (x^/ndf « 0.992). The total x^/ndf for ah three processes is 
379.893/(301 - 42) ^ 1.46. 

The background parameters are an = 0.183, aia- = 0.0252, ai^ = 0.0155, 
bu = 0, bia = -0.0089, biy = 0.04336, a2i = -0.6973, a2a = -1.427, a2v = 
-5.935, 621 = 0.0447, = 0, = 7.044, 631 = 0.6346, ha = 0.3336, 
62^ = 0;sa = 1.638 GeV^ = 2.084 GeVl 
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On figures 1-5, we demonstrate energy dependences of phase shifts and 
moduh of the matrix elements of processes tttt — ^ tttt, KK^ rjrj compared 
with the experimental data. 
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Figure 1: The phase shift of the TTTT-scattering 5"- wave ampUtude. 
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Figure 2: The module of the TTTT-scattering S-wscve matrix element. 
Let us indicate in Table 1 the obtained pole clusters for resonances 
on the complex energy plane -s/s, poles on sheets IV, VI, VIII and V, 
corresponding to the /o(1500), are of the 2nd and 3rd order, respectively 
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7T + TT ^ K + K 
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Figure 3: The phase shift of the tttt KK S-w&ve matrix element, 
(this is an approximation). 

Table 1: Pole clusters for the /o-resonances. 



Sheet 


II 


III 


IV 


V 


VI 


VII 


VIII 


/o(600) 


r. 


678±14 
608±22 


688±16 

608±9 






628±17 

608±28 


618±15 

608±26 




/o(980) 


E^ 

r. 


1016±5 
32±8 


986±18 
59±16 












/o(1370) 


E^ 

r. 








1400±21 
89±13 


1400±20 
71±15 


1400±20 
45±6 


1400±20 
27±9 


/o(1500) 


E^ 

r. 


1505±22 
360±23 


1480±30 
140±21 


1505±20 
240±30 


1500±20 
139±21 


1493±27 
194±27 


1488±25 
88±15 


1505±20 
360±30 


/o(1710) 


E^ 




1704±18 
95±14 


1704±21 
105±17 


1704±32 
325±26 


1704±30 
325±45 







Note a surprising result obtained for the /o(980) state. It turns out that 
this state hes slightly above the KK threshold and is described by a pole on 
sheet II and by a shifted pole on sheet III under the ijij threshold without 
an accompaniment of the corresponding poles on sheets VI and VII, as it 
was expected for standard clusters. This corresponds to the description of 



7T + 71 ^ K + K 
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Vs" [MeV] 

Figure 4: The module of the tttt — > KK S-wme matrix element. 



the T]r] bound state. 

For now, we did not calculate coupling constants in the 3-channel ap- 
proach. Therefore, for subsequent conclusions, let us mention the results 
for coupling constants from our previous 2-channel analysis (Table 2) [15]: 
gir is the coupling constant of resonance "r" with the TTTr-system; g2r, with 
KK. We see that the /o(980) and the /o(1370) are coupled essentially 



Table 2: Coupling constants of the /o-resonances from the 2-channel analysis. 





/o(665) 


/o(980) 


/o(1370) 


/o(1500) 


9ir. GeV 


0.652 ± 0.065 


0.167 ±0.05 


0.116 ±0.03 


0.657 ±0.113 


92r, GeV 


0.724 ±0.1 


0.445 ±0.031 


0.99 ±0.05 


0.666 ±0.15 



more strongly to the KK system than to the tttt one, i.e., they have a dom- 
inant ss component. The /o(1500) has the approximately equal coupling 
constants with the tttt and KK, which apparently could point to its dom- 
inant glueball component. In the 2-channel case, /o(1710) is represented 
by the cluster corresponding to a state with the dominant ss component. 
Our 3-channel conclusions on the basis of resonance cluster types gen- 
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71 + 71 ^ ?7 + ?7 
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Figure 5: The squared module of the tttt rjr] S-wave matrix element. 

erally confirm the ones drawn in tlie 2-cliannel analysis (besides the above 
surprising conclusion about the /o(980) nature). 

Masses and widths of these states should be calculated from the pole 
positions. If to take the resonance part of amplitude as 



we obtain for masses and total widths the following values (in MeV) : 

for /o(600), 868 and 1212; 

for /o(980), 1015.5 and 64; 

for /o(1370), 1407.5 and 344; 

for /o(1500), 1546 and 716; 

for/o(1710), 1709.6 and 276. 



(10) 



4 Lower scalar nonets 

It is known that an assignment of the scalar mesonic states to quark- 
model configurations is problematic up to now, although there is a number 
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of interesting conjectures [16]- [21]. Let us also (on the basis of obtained 
results) propose a following assignment of scalar mesons below 1.9 GeV to 
lower nonets. First of all, we exclude from this consideration the /o(980) 
as the ijij bound state. Then we propose to include to the lowest nonet the 
isovector ao(980), the isodoublet KQ^dOb) (or /^(800)), and two isoscalars 
/o(600) and /o(1370) as mixtures of the eighth component of octet and the 
SU(3) singlet. Note that we consider the Kq{905) (or k) which one has 
observed at analysing the K — tt scattering [22, 23], extracted from reaction 
K~p — ^ K~7r^n, and at studying the decay — ^ K'tt^tt^ [24]. Then 
the Gell-Mann-Okubo formula 

3m/g = ArriK* - rua^ (11) 

gives TTLf^ = 0.88 GeV. Our result for the cr-meson mass is rria- ~ 0.868±0.02 
GeV (if m« = 0.8, m/, « 0.73). 

In the relation for masses of nonet 

ma + m/Q(1370) = 2m7^*(905)), (12) 

the left-hand side is about 26 % bigger than the right-hand one if to take 
our mass values. 

The next nonet (maybe, of radial excitations) could be formed of the 
isovector ao(1450), the isodoublet A'o*(1430 - 1450), and of the /o(1500) 
and /o(1710) as mixtures of the eighth component of octet and the SU(3) 
singlet, the /o(1500) being mixed with a glueball which is dominant in this 
state. From the Gell-Mann-Okubo formula we obtain nif^ 1.45 GeV. In 
second formula 

"^/o(1500) + "^/o(1710) = 2m;i'*(i45o), (13) 

the left-hand side is about 12 % bigger than the right-hand one. 
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Though the Gell-Mann-Okubo formula is fulfiUed for both nonets rather 
satisfactorily, however, the breaking of 2nd relation (especially for the low- 
est nonet) tells us that the cr — /o(1370) and /o(1500) — /o(1710) systems get 
additional contributions absent in the i^o(905) and A'o(1450), respectively. 

Note that the 3-channel analysis indicates on a non-simple picture of 
mixing the states /o(1370) and /o(1710) with the wide /o(1500) and /o(600). 
The /o(1370) is coupled more strongly to the 7777 channel than to the tttt 
and KK ones; the /o(1710), to the KK than tttt and 7777 ones, whereas if 
these states were the pure ss ones, they would be described by clusters of 
type (e) , and their coupling constants with the KK and 7777 channels would 
be congruent numbers for each state. 



5 Model analysis of isoscalar-tensor sector 



We analyze the isoscalar D-waves of the processes tttt — ^ tttt, KK, rir] in 
the 4-channel approach with the explicit account of the channel (27r) (27r) 
(i=4), too. The Jost matrix determinant d{ki, k2, k^, k^) is taken as 



d — d^d, 



res- 



(14) 



The 4-channel Breit-Wigner form for the resonance part of the (i-function 



is taken in the form {prj = 2ki/^M^ — 4mj): 



dres{s) = n 



M,' - S - i E PljRrjfrj 



(15) 



where f^j/Mr is the partial width, Rrj is a Blatt-Weisskopf barrier factor 
[8] 



9 + K^s - 4mj Trjf + - ^rjY 
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with radii of 0.955 Fermi for all resonances in all channels, except for 
/2(1270), /2(1525) and /2(1950) for which they are: for /2(1270), 1.496, 
0.704 and 0.604 Fermi respectively in channels tttt, KK and rjr] , for 
/2(1525), 0.576 and 0.584 Fermi in channels KK and 7777, and for /2(1950), 
0.178 Fermi in channel KK. 

The background is parameterized by 



ds = exp 



/ 2/Cr 



n=l V V S 



(17) 



To take into account an influence of opening channels in the range of ~ 
1.5 GeV (7777', pp, iJij), ai and 6„ (n = 1,2,3) are taken in the form: 

ai = an H —ai29(s - 4mi) H ^^lo^'fs - Sy)), (18) 

s s 

s — s 

bn = ^n-\ --fnO{s-Sy), (19) 

s 

where Sy ^ 2.274 GeV^ is the combined threshold of channels rjrj', pp, uu. 

The data for the tttt scattering are taken from an energy-independent 
analysis by B. Hyams et al. [11]. The data for the processes tttt — ^ KK, rjr] 
are taken from works [25]. We obtained ten resonances /2(1270), /2(1430), 
/^(1525), /2(1580), /2(1730), /2(1810), /2(1950), /2(2000), /2(2240) and 
/2(2410). We will not discuss the last resonance because there are not 
practically data. We need it for satisfying unitarity. 

On figures 6-9, we demonstrate results from our fitting to data (in the 
Argand plot units). 



We obtain a reasonable description (the total x^/ndf = 162.577/(168 — 
66) ~ 1.59) with the values of parameters of /2-resonances shown in Table 
3. 
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7T + 71 ^ 71 + 7T 




750 1000 1250 1500 1750 2000 2250 

Vs [MeV] 

Figure 6: The phase shift of the TTTr-scattering D-wave amphtude. 



For the background we find: 

ail = -0.0785, ai2 = 0.0345, = -0.2342, Pi = -0.06835, 

71 = -0.04165, /?2 = -0.981, 72 = 0.736, /Ss = -0.5309, 73 = 0.8223. 

6 Discussion and conclusions 

• In combined 3-channel model-independent analysis of data on pro- 
cesses TTTT — ^ TrTTjKKjrjr] in the channel with I'^J^'^ = 0"^0"^"^, an 
additional confirmation of the cr-meson with mass 0.868 GeV is ob- 
tained. This mass value rather accords with prediction [rrifj « nip) on 
the basis of mended symmetry by Weinberg [26]. In works [27]-[32] 
evidences of the existences of the cr-meson have been given, too. 

• The /o(1370) and /o(1710) have the dominant ss component. Con- 
clusion about the /o(1370) quite well agrees with the one of work of 
Crystal Barrel Collaboration [33] where the /o(1370) is identified as 7777 
resonance in the Tr^rji] final state of the pp annihilation at rest. Con- 
clusion about the /o(1710) is quite consistent with the experimental 
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7T + 7T ^ 7T + 7T 
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Figure 7: The module of the TTTr-scattering D-wave matrix element. 

facts that this state is observed in 77 — ^ KsKg [34] and not observed 
in 77 TT^TT^ [35]. 

• Indication for /o(980) to be the rji] bound state is obtained. From 
point of view of quark structure, this is the 4-quark state. Maybe, 
this is consistent somehow with arguments of Refs. [36, 37] in favour 
of the 4-quark nature of /o(980). 

Remembering a dispute [38, 39] whether the /o(980) is narrow or not, 
we agree rather with the former. Of course, it is necessary to make 
analysis of other relevant processes, first of all, J/il) and (/) decays. 

• As to the /o(1500), we suppose that it is practically the eighth compo- 
nent of octet mixed with a glueball being dominant in this state. Its 
biggest width among enclosing states tells also in behalf of its glueball 
nature [18]. 

• An assignment of the scalar mesons below 1.9 GeV to lower nonets is 
proposed. Note that this assignment moves a number of questions and 
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7T + 71 ^ K + K 
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Figure 8: The squared module of the tttt — > KK D-wave matrix element. 



does not put the new ones. Now an adequate mixing scheme should 
be found. 

• We do not obtain /2(1640), /2(1910), /2(2150), /2(2010), however, we 
see /2(1450) and /2(1730) which are related to the statistically- valued 
experimental points. 

• Usually one assigns to the first tensor nonet the states /2(1270) and 
/2(1525). To the second nonet, one could assign /2(1601) and /2(1767) 
though for now the isodoublet member is not discovered. If one takes 
for the isovector of this octet the state a2(1730) and if the /2(1601) 
is almost its eighth component, then, on the basis of the Gell-Mann- 
Okubo formula, we would expect this isodoublet mass at about 1.635 
GeV. Then the relation for masses of nonet would be well fulfilled. 
Note that in the Particle Data Group issue [6] is indicated an exper- 
iment [40] in which one has observed the strange isodoublet with yet 
indefinite remaining quantum numbers and with mass 1.629 ± 0.007 
GeV. 
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Figure 9: The squared module of tttt — > 7777 D-wave matrix element. 

• The states /2(1963) and /2(2207) together with the isodoublet i^:2*(1980) 
could be put into the third nonet. Then in the relation for masses of 
nonet 

-^/2(1963) + -^/2(2207) = 2Mi;j:*(i98o), (20) 

the left-hand side is only 5.3 % bigger than the right-hand one. If one 
consider /2(1963) as the eighth component of octet, then the Gell- 
Mann-Okubo formula 

Ma, = 4M^*(l980) - 3Mj^(i963) (21) 

gives = 2.031 GeV. This value practically coincides with the one 
(2.03 GeV) for a2-nieson obtained on the basis of the recent data [41]. 
This state is interpreted [7] as a second radial excitation of the 1~2++- 
state on the basis of consideration of the a2 trajectory on the (n, M^) 
plane (n is the radial quantum number of the qq state) . 

• As to 72(2017), the ratio of the tttt and rjr] partial widths is in the 
limits obtained in Ref.[7] for the tensor glueball on the basis of the 



18 



Table 3: The /2-resonance parameters (all in the MeV units). 



State 


M 


frl 


fr2 


fr3 


/r4 




/2(1270) 


1275.1±1.8 


470.9±5.4 


201.5±11.4 


89.5±4.76 


22.6±4.6 


>212 


/2(1430) 


1450.8±18.7 


128.3±45.9 


562.3±142 


32.7d=18.4 


8.2±65 


>230 


/^(1525) 


1535±8.6 


28.6±8.3 


253.8±78 


92.7±11.5 


41.4±160 


>76 


/2(1565) 


1601.4±27.5 


75.5±19.4 


315±48.6 


388.9±27.7 


127±199 


>170 


/2(1730) 


1724.4±5.7 


78.8±43 


289.5±62.4 


460.3±54.6 


107.6±76.7 


>181 


/2(1810) 


1766.5±15.3 


129.5±14.4 


259±30.7 


469.7±22.5 


90.3±90 


>177 


/2(1950) 


1962.8±29.3 


132.6±22.4 


333±61.3 


319±42.6 


65.4±94 


>119 


/2(2000) 


2017±21.6 


143.5±23.3 


614±92.6 


58.8±24 


450.4±221 


>299 


/2(2240) 


2207±44.8 


136.4±32.2 


551±149 


375±114 


166.8±104 


>222 


/2(2410) 


2429±31.6 


177±47.2 


411±196.9 


4.5±70.8 


460.8±209 


>169 



1/N-expansion rules. However, the KK partial width is too large for 
the glueball. This question requires an additional investigation. 

• Finally we have /2(1450) and /2(1730) with the rather unusual prop- 
erties. These are non-qq states and non-glueball. Since one predicts 
that masses of the lightest qqg hybrids are bigger than the ones of 
lightest glueballs, maybe, these states are the 4-quark ones. 
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Program. M.N. were supported in part by the Slovak Scientific Grant 
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